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Abstract 

The precise measurements of the B® oscillation frequency and the limit on the B® one 
as well as the determination of the Cabibbo-Kobayashi-Maskawa matrix element |V u t>| 
improve the constraints on the other elements of this matrix. 

A fit to the experimental data and the theory calculations leads to the determination 
of the vertex of the unitarity triangle as: 

„ _ n i cn+0.094 n 001 +0.061 

p — U.lDU_ o7o V — u -"J oi -0.058- 

The values of its angles, in their customary definition in terms of sines for a and /?, are 
found to be: 

sin2a = 0.06ig;l sin 2/5 = 0.75 ± 0.09 7 = 67±^°. 

Indirect information on non-perturbative QCD parameters, on the presence of a CP 
violating complex phase in the CKM matrix and on the B® oscillation frequency are also 
extracted. 

To be submitted to Physics Letters B 
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Introduction 



The Standard Model P of the electroweak interactions predicts a mixing of the quark mass 
eigenstates with the weak interaction ones. This mixing is described by the Cabibbo-Kobayashi- 
Maskawa jSj (CKM) matrix. Four real parameters describe this 3x3 unitary matrix j^j: 

/V ud V us V ub \ / 1-f A AX 3 (p- iV ) 

V cd V cs V cb ] ~ I -A 1-f AX 2 |+0(A 4 ). (1) 
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CKM 



\v td v ts v tb y v^A 3 (i-p-^) —ax 



As A, p and r] are of order unity, and A is chosen as the sine of the Cabibbo angle, this 
parametrisation shows immediately the hierarchy of the couplings of the quarks in the charged 
current part of the Standard Model Lagrangian. Moreover in this parametrisation the parameter 
7] is the complex phase of the matrix and is thus directly related to the known violation of the 
CP symmetry produced by the weak interactions. The measurement of the parameters of the 
CKM matrix is thus of fundamental importance for both the precision description of the weak 
interaction of quarks and the investigation of the mechanism of CP violation. 

The parameters A and A are known with an accuracy of a few percent and the determination 
of p and rj is the subject of this letter. A large number of physical processes can be parametrised 
in terms of the values of the elements of the CKM matrix, together with other parameters 
of theoretical and experimental origin. Four of them show good sensitivity for the indirect 
determination of p and rj and are discussed in what follows. A fit based on this information 
has been performed, as suggested in [5], and its results are presented below. 

As it is well known the measurement of p and r] is equivalent to the determination of the 
only unknown vertex and the angles of a triangle in the p — i] plane whose other two vertices 
are in (0,0) and (1,0). Figure 1 shows this triangle, called the unitarity triangle. 

Constraints 

The value of the sine of the Cabibbo angle is known with a good accuracy as: 

A = 0.2196 ±0.0023. 

The parameter A depends on A and on the CKM matrix element |V c b|. Using the value 

|Vcb| = (39.5 ±1.7) x 10~ 3 , 



A = = 0.819 ±0.035. 



it can be extracted: 

A 2 

The four processes most sensitive to the value of the CKM parameters p and r\ are described 
in the following, along with their experimental knowledge and theoretical dependences. 



CP Violation for Neutral Kaons 

The mass eigenstates of the neutral kaons can be written as \K$) = p\K°) ± q\K°) and \Kl) = 
p\K°) — q\K°). The relation p ^ q implies the violation of CP that, in the Wu-Yang phase 
convention is described by the parameter ex defined as: 

p _ l + e K 
q 1 - e K ' 
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The precise measurements of the Kg — > 7r + 7r and Kl — ► vr + 7r decay rates imply [Hj: 

\e K \ = (2.280 ±0.019) x 10~ 3 . 
The relation of \ek\ to the CKM matrix parameters is [HIE]: 
G 2 F flm K m 2 w 



\ e K\ 



-B K (A 2 \ G ri) [y c ( Vct f 3 (y c , y t ) - Vcc ) + Vuytf2(y t )A 2 \ 4 (1 - p) ]. (2) 



6y/2TT 2 Am K 

The functions / 3 and / 2 of the variables y t = m 2 /m 2 v and ?/ c = m 2 /m^ are given by 
, / x 1 9 3 3x 2 In a; y 3y / ylnj/\ 

/2(x) = i + 4(i^-(r^-2(r^F' /3(a; ' y) = ln x-4(r^) ( 1 + i^J- (3) 

From the value of the mass of the top quark reported by the CDF and DO collaborations [Sj, 
173.8 ± 5.2 GeV, and the scaling proposed in [H] one obtains: 

fn t (m t ) = 166.8 ± 5.3 GeV, 
while the mass of the charm quark is [5] : 

nr c (m c ) = 1.25 ±0.15 GeV. 

The current values of the parameters that include the calculated QCD corrections are |9*lll0j: 

r] cc = 1.38 ± 0.53, rjtt = 0.574 ±0.004, r] ct = 0.47 ± 0.04. 

The largest theoretical uncertainty which affects this constraint is that on the "bag" parameter 
Bk, that reflects non-perturbative QCD contributions to the process. Using the value of the 
JLQCD collaboration ^T], -B_r-(2 GeV) = 0.628 ± 0.042, with a calculation similar to that 
reported in tQj the value used in the following can be derived as: 

B K = 0.87 ±0.14. 

The other physical constants of the formula are reported in Table 1. This constraint has 
the shape of an hyperbola in the p — i] plane. 



Oscillations of Mesons 

Neutral mesons containing a b quark show a behaviour similar to neutral kaons. The heavy 
and light mass eigenstates, B L and B H respectively, are different from the CP eigenstates B% 
and B° d : 

\B L ) = p\B° d ) + q\Bl) \B H ) = p\B° d ) - q\B°). 

In the neutral B system the mass difference Am^ = mg fl — m# L is the key feature of the 
physics while the lifetime difference dominates the effects in the neutral kaon system. This 
mass difference can be measured by means of the study of the oscillations of one CP eigenstate 
into the other. The high precision world average is |13j : 

Am d = 0.471 ±0.016 ps" 1 . 
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The relation of Ama with the CKM parameters, making use of the Standard Model description 
of the box diagrams that give rise to the mixing and the parametrisation (0) of the CKM 
matrix, reads: 

Am d = ^m 2 w m B (f^^B^ 2 VB y t f 2 (y t )AX 6 [(1 - pf + r, 2 ] . (4) 

The function f 2 is given by ©, the value of the calculated QCD correction r\ B is PITU]: 

T) B = 0.55 ±0.01, 



and the equivalent of the B K parameter for the kaon system, f Bd JB Bd , is taken as 



f Bd ^B Bd = 0.201 ±0.042 GeV. 

The measurement of Amj constrains the vertex of the unitarity triangle to a circle in the 
p — r] plane, centred in (1, 0). 

Oscillations of Mesons 

B® mesons are believed to undergo a mixing analogous to the B® ones. Their larger mass 
difference Am s is responsible for oscillations that are faster than the B® ones, and have thus still 
eluded direct observation. A lower limit has been set by the LEP B oscillation working group 
combining the results of the searches performed by the LEP experiments with a contribution 
from the SLD and CDF collaborations, as j!5j : 

Am s > 12.4 ps- 1 (95%C.L.). 

The expression for Am s in the Standard Model is similar to that for Am^. From the ratio 
of these two expressions the value of Am s can be written as: 

Am, = AmA^Z 2 - \- 2 -, (5) 

\ 2 m Bd (l- p ) 2 + r? 2 

where all the theoretical uncertainties are included in the quantity £, known as [Tlj : 



fs d \ B B 
£ = V ' = 1.14 ±0.08. 

JB 3 V B Bs 

This experimental lower limit excludes the values of the vertex of the unitarity triangle 
outside a circle in the p — t] plane with centre in (1,0). 

Charmless Semileptonic b Decays 

The three constraints described above are all affected by a large theoretical uncertainty on some 



of the parameters that enter their expression, namely Bk, f Bd y B Bd and £. A determination 
of either |V u b| or the ratio |V u b|/|V c b| allows a more sensitive constraint not relaying on any 
non-perturbative QCD calculation. It follows from the CKM matrix parametrisation of (fljl 
that: 

|V ub |/|V cb | = \^p 2 + v 2 - (6) 



4 



The CLEO collaboration has measured this ratio by means of the endpoint of inclusive [T6] 
charmless semileptonic B decays as: |V u b|/|V c b| = 0.08 ± 0.02. The ALEPH and L3 collabo- 
rations have recently measured at LEP the inclusive charmless semileptonic branching fraction 
of beauty hadrons, Br (b — ► X u £i/), from which the value of |V u b| can be extracted as: 

/Br(6 -»• XJv) /l.6ps 
|V ub | = 0.00458 x y 1 ooQ2 M ; x J— -Z- ± 4% thcory . (7) 

The experimental results are : 

ALEPH [18 : Bi(b -> XJz/) = (1.73 ± 0.55 ± 0.55) x 1(T 3 
L3 [H]: Br(6 -> XJz/) = (3.3 ± 1.0 ± 1.7) x lO" 3 , 

where the first uncertainty is statistical and the second systematic, with the average: 

Br(6 -> XJz/) = (1.85 ± 0.52 ± 0.59) x 1(T 3 , 

with the same meaning of the uncertainties. This value makes it possible to determine |V u b| at 
LEP by means of the formula as: 

|V u b| = (4.5l°;^i±0.2)xl(r 3 . 

The first uncertainty is statistical, the second systematic and the third theoretical. The value 
r B = (1.554 ± 0.013)ps [20] has been used. Using the quoted value of |V c b| the combination 
with the CLEO measurement yields: 

|V ub |/|V cb | = 0.093 ±0.016. 

The uncertainty on this important constraint is thus significantly reduced by the inclusion of 
the recent LEP measurements. A further reduction to 0.015 could be achieved by the inclusion 
of the DELPHI collaboration preliminary measurement of this quantity |21j . 

This constraint gives a circle in the p — r\ plane with centre in (0,0), shown in Figure 2 
together with all the other constraints described above. 

Determination of p and 77 

The p and rj parameters can be determined from a fit to the experimental values of all the 
constraints described above. The experimental and theoretical quantities that appear in the 
formulae describing the constraints have been fixed to their central values if their errors were 
reasonably small, and are reported in the left half of Table 1. The quantities affected by a 
larger error have been used as additional parameters of the fit, but including a constraint on 
their value. This procedure has been implemented making use of the MINUIT package (22j to 
minimise the following expression: 
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Fixed in the fit 


Varied in the fit 


\ 

A 


n 01 net _i_ n nnoo 
= U.zlyb ± U.UUzo 


A 


n 010 _i_n nop; 
= U.oly ± U.UoD 




= (i.ioooy ± u.uuuuij x iu oev 


r]ct 


n a 1 j- n n/i 
= U.4 / ± U.U4 


JK 


n 1 c^Qs -i- n nni ^ p v 

— u.ioyo it u.uuio vjrev 


Vcc 


— 1.00 it u.oo 


Amj/- 


= fO 5304 ± 0014) x 10 _2 ds _1 


Tn~( rn,„ ) 


= 1 25 ± 15 GeV 


TflK 
J\ 


= 0.497672 ± 0.000031 GeV 


rfu(rrit) 


= 166.8 ±5.3 GeV 








n 9D1 4- D 0.49 HpV 


m Bd 


c 0700 -i- n nni s n„u 
= o.z ii)z ± u.uuio Lrev 




n Q7 1 n n 
= U.o 1 ± U. 14 


m B s 


= o.ooyz ± U.UUzU Lrev 




1 i-i in ns 
= 1.14 ± U.Uo 


m B 


— o.zyu m u.uuz vjev 




(9 9sn 4- n m q\ v 1 n -3 


Vb 


= 0.55 ±0.01 




= 0.471 ±0.016 ps- 1 


Vtt 


= 0.574 ±0.004 


|V ub |/|V cb | 


= 0.093 ±0.016 



Table 1: Physical constants and parameters of the fit. The values whose origin is not discussed 
in the text are from [5]. 



(W| - |e^|) [Am d -Am d 



+- - 2 — + A - 2 ± X (A {Am.) , a A (Am,)) . 

a \e K \ a Am d 

The symbols with a hat represent the reference values measured or calculated for a given 
physical quantity, as listed in Table 1, while the corresponding a are their errors. The param- 



eters of the fit are p, r], A, m c , m t , Bk, rj ct , i] cc , fB d \JBB d and £, that are used to calculate the 
values of |e#|, Am d , Am s and by means of the formulae (0), Q, © and ©• 

As no measurements of Am s are available a further contribution to the \ 2 analogous to the 
previous ones can not be calculated. The following approximation has been used to extract a 
contribution from the Confidence Levels of the Am s exclusion. The results of the search for B° s 
oscillations have been presented and combined ^H] in terms of the oscillation amplitude A [23J, 
a parameter that is zero in the absence of signal and compatible with one if an oscillation signal 
is observed, as in: 



P 



B° s -> (B°,B°)] = ^-e~ t/Ts (1 ±^cosAm s 
J 2t k 



The experimental results are reported in terms of A(Am s ) and 04 (Am,), which leads to the 
quoted 95% Confidence Level limit as the value of Am s for which the area above one of the 
Gaussian distribution with mean A (Am,) and variance a\ (Am,) equals the 5% of the total 
area. As noted in |23J the full set of combined A (Am,) and oj± (Am,) measurements indeed 
contains more information than this limit and it is used in this procedure, with a different 
statistical approach. The value of Am s can be calculated for each value taken by the fit 
parameters p, 77 and £ by means of formula (JSJ), together with the value of its corresponding 
Confidence Level obtained as described above. The value \ 2 (-4 (Am,) , 04 (Am,)) of a \ 2 
distribution with one degree of freedom corresponding to this Confidence Level can then be 
calculated and added to the total \ 2 °f the fit. 
The results of the fit are the following: 

„_ni cn+0-094 _ n 001 +0.061 

p - 0.160_ 0070 7] - 0.3b4„ 0058 . 
The 95% Confidence Level regions for p and r\ are: 

-0.02 < p < 0.35 0.27 < 77 < 0.50 (95%C.L.). 
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Figure 3a shows the allowed confidence regions in the p — rj plane, together with the favoured 
unitarity triangle, that is also shown superimposed on the constraints of Figure 2. 

From these results it is possible to determine also the value of the angles of the unitarity 
triangle. The angles a and (3 are reported in terms of the functions sin 2a and sin 2(3 as will be 
measured at the next B-factories. The numerical values obtained from the fit are: 

sin 2a = 0.06±g;^ sin 2(3 = 0.75 ± 0.09 7 = 67t^°. 

In terms of 95% Confidence Level regions these last results can be expressed as: 

-0.71 < sin 2a < 0.70 0.56 < sin 2(3 < 0.94 44° < 7 < 93° (95%C.L.) 



Interpretations 

The fit procedure described above can also be used to extract information on the theory param- 
eters that enter the fit with a large uncertainty and at the same time, perform an estimation of 
p and 77 independent of them. This can be achieved by removing from the fit the constraint on 



the parameter. The two parameters B K and fB d ^jB Bd are those affected by the largest theory 
uncertainty. By applying this method to the parameter Bk, the fit yields: 

p = 0.156lffi ? = 0.393j$:K B K = 0.W%%. 

The value of Bk favoured by the fit has an error larger than that on the estimated input 
parameter and thus can not help in restricting its range of allowed values. The same procedure 



with fB d \ Bs d as a free parameter leads to the results: 



p = 0.186±8:8ii 77 = 0.379^ = - 222 ±oSu GeV 



the value of fB d \jB Bd comes out to be well in agreement with the predicted one with a smaller 
uncertainty. The same procedure applied to B K and fB d \B Bd simultaneously gives: 



p = 0.170l£S£ 77 = 0.3901^3 B K = 0.82%£ f B ^B Bd = 0.217^ GeV. 

The Am s constraint has a big impact on the p uncertainty as can be observed by removing 
it from the fit, what gives: 

a - 022+ - 145 77-O 4S4+ - 062 

Figure 3b shows the experimentally favoured regions in the p — 77 plane for this fit together with 
the lower limit and expected sensitivity (Am s = 13.8 ps _1 [T^j) of the current experiments to 
B° s oscillations. The confidence regions for Am s can be extracted from this fit as: 

Am s = 11.3t|£ ps -1 

5.7PS" 1 < Am, < 17.8 ps" 1 (95%C.L.), 

The LEP measurements have greatly improved the constraints on the CKM matrix. Another 
fit has been performed removing the Am s constraint, derived mainly from the LEP limits, and 
excluding the LEP measurement from the averages of the other input quantities; that is using: 

|V ub |/|V cb | = 0.08 ±0.02 

|V cb | = (50±5)xl0~ 3 
Am d = (0.500 ± 0.030) ps" 1 . 
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The first value is that quoted above from the CLEO collaboration jTHj, the second follows 
from and the last has been estimated from the current published and preliminary results 
from the CDF and SLD collaborations. This fit, as shown in Figure 3c, yields: 

p = 0.012l°il r) = 0.3831K and 

sin 2a = 0.63±g;|5 sin 2(3 = 0.67j#M 7 = 88±|| °. 
Some of the errors are reduced by as much as a factor three by the inclusion of the LEP data. 

A real CKM matrix ? 

To date the only experimental evidence for the violation of CP in the CKM matrix, namely 
its complex phase described by a value of r\ different from zero, comes from the neutral kaon 
system. As different models have been proposed to explain that effect, it is of interest to remove 
from the fit the constraint related to this process and then investigate the compatibility of r\ 
with zero |2E]- This procedure yields the following results, graphically displayed in Figure 3d: 

n _ n 1 c; fi +0.096 q oq 4 +0.072 

P — u - iou -0.090 7 — u - oy ^-0.080- 

The value of rj is not compatible with zero at the 95% and 99% of Confidence Levels either: 

-0.025 < p < 0.358 (95%C.L.) -0.069 < p < 0.411 (99%C.L.) 
0.224 < r] < 0.531 (95%C.L.) 0.157 < r? < 0.574 (99%C.L.). 

If the CKM matrix is assumed to be real, as recently proposed for instance in [2Zj, all the 
circular constraints reduce to linear intervals on the p axis, onto which the unitarity triangle will 
then be projected. This hypothesis can be checked removing again the neutral kaon constraints 
from the fit and modifying the formulae (J1J, (J5J) and © imposing r\ equal to zero. The result 



of this fit, whose parameters are reduced to p, A, m t , fB d \jBB d and £, is: 

n - n QOT +0-053 
P — u - ozi -0.056- 

The value of the x 2 function at the minimum is 6.7, leading to the conclusion that a CKM 
matrix real by construction can fit the data. 



Conclusions 

The combination of the precise measurements of Am^, the updated limits on Am s and the 
determination of |V u b| helps in constraining the CKM matrix elements. 

From a simultaneous fit to all the available data and theory parameters the vertex of the 
unitarity triangle is determined as: 

„_ni cn+0-094 _ n oo-i +0.061 

p - 0.160_ 070 7] - 0.381_ 058 . 
yielding the following values for its angles: 

sin 2a = 0.06^2 sin 2(3 = 0.75 ± 0.09 7 = 67t^°. 
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The accuracy on sin 2(3 from this indirect analysis is already at the same level as that ex- 
pected to be achieved with the direct measurement at the B-factories due to become operational 
in the next future. These limits greatly benefit from the inclusion of LEP data. 

The fit suggests the value of the non-perturbative QCD parameter f Bd \jB Bd as: 

^ V / ^ = - 222 -o:on GeV. 

The parameter 7] related to the complex phase of the matrix and thus to the CP violation 
is found to be different from zero at more than the 99% Confidence Level, even removing from 
the fit the constraints arising from CP violation in the neutral kaon system. Nonetheless the 
hypothesis of a real matrix can still fit the data without this constraint. 

The fit also indicates the Am s variation range as: 

Am s = (11.3+|°) ps" 1 

5.7 ps" 1 < Am s < 17.8 ps" 1 (95%C.L.). 

These results improve those of similar previous analyses [2UEB] an d agree with another 
one based on a different approach [TK] . 
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(0,0) 




(1.0) p 



Figure 1: The unitarity triangle. 




Figure 2: The current constraints and the favoured unitarity triangle. The constraint coming 
from B° oscillations is a limit at 95% of Confidence Level, while the others represent a ±lcx 
variation of the experimental and theoretical parameters entering the formulae in the text. 
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p p 



Figure 3: The favoured unitarity triangles and the confidence regions for their vertices in the 
following assumptions: a) the fit using all data described in the text, b) the constraint from 
the B° oscillations is not applied, c) the LEP measurements are excluded from the fit, d) no 
constraints from the neutral kaon system are applied. The band in b) displays the values of p 
and rj corresponding to a value of Am s between the current lower limit and expected sensitivity. 
The Am s limit and the central values of the constraints are shown in a), c) and d). 
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